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We study, within a framework of the classical fields approximation, the density correlations of a 
weakly interacting expanding Bose gas for the whole range of temperatures across the Bose-Einstein 
condensation threshold. We focus on elongated quasi-one-dimensional systems where there is a 
huge discrepancy between the existing theory and experimental results (A. Perrin et ah. Nature 
Phys. 8, 195 (2012)). We find that the density correlation function is not reduced for temperatures 
below the critical one as it is predicted for the ideal gas or for a weakly interacting system within 
the Bogoliubov approximation. This behavior of the density correlations agrees with the above 
mentioned experiment with the elongated system. Although the system was much larger then studied 
here we believe that the behavior of the density correlation function found there is quite generic. Our 
theoretical studies indicate also large density fluctuations in the trap in the quasicondensate regime 
where only phase fluctuations were expected. We argue that the enhanced density fluctuations 
can originate in the presence of interactions in the system, or more precisely in the existence of 
spontaneous dark solitons in the elongated gas at thermal equilibrium. 


Correlations are the essence of any many body sys¬ 
tems. In particular, the quantum features of the system 
are manifested by some unusual correlations. The first 
order coherence is the basic criterion used as the defi¬ 
nition of a Bose-Einstein condensate. And indeed, very 
soon after observation of trapped atomic condensates the 
first order coherence of such systems, manifesting itself in 
the ability to produce the interference fringes in a two-slit 
experiment, had been proven [1] experimentally. Higher 
order correlations leading to atom bunching were estab¬ 
lished from collisions and three-body losses [2]. Although 
the Glauber coherence theory introduced to characterize 
correlations of quantum electromagnetic field is well es¬ 
tablished now, the issue of coherence of a matter field is 
still under intensive investigation. 

A Bose-Einstein condensate is a matter wave analogue 
of a coherent light. A very important question is how far 
this analogy can be pursuit. The first difference is that 
atoms exist in a Fock states only: the states which are the 
eigenstates of the particle number operator. As a conse¬ 
quence, the coherent state of a matter field, understood 
as the exact analogue of the coherent electromagnetic 
field, does not exist. On the other hand, the atomic field 
can exhibit a higher order coherence too. The coherence 
ought to be understood here as the ability to produce 
the interference patterns not only in a single-particle de¬ 
tection schemes, but also in a simultaneous detection of 
larger number of atoms. Although the phase of a num¬ 
ber state can be arbitrary, the Fock states can interfere 
in a single realization of the system [3] when no averaging 
over global phase is performed. Therefore higher order 
coherence of atomic condensates, as can be observed in 
a single realization of the system, is an interesting issue 
and in fact this is the property of the system which char¬ 
acterizes a genuine Bose-Einstein condensate [4]. The 
correlation functions of a Bose-Einstein condensate - a 


coherent atomic wave - can be contrasted with the corre¬ 
lations of fermionic atoms and/or thermal nondegenerate 
clouds. 

The famous experiment by R. Hanbury Brown and 
R.Q. Twiss [5] exploring a new type of interferometer 
to measure the angular diameter of radio stars initiated 
a modern theory of coherence. The Hanbury Brown and 
Twiss effect is rooted in properties of the density-density 
correlations or equivalently, the second order correlation 
function. The bunching effect for thermal photons origi¬ 
nating from many particle interference enhancing proba¬ 
bility of many photons detection can be contrasted to 
the absence of such an effect for coherent light when 
only a single mode is populated. The atom counting 
experiments with metastable bosonic or fermionic he¬ 
lium atoms, exploring thermal density-density correla¬ 
tions, show bunching of nondegenerate bosons and an¬ 
tibunching effects for fermions [6]. Similar results were 
obtained in experiments with nondegenerate gas at opti¬ 
cal lattices [7] and atomic lasers [8]. Direct measurement 
of the second and the third order correlation function of 
Bose-Einstein condensates is an important proof of their 
high coherence [9]. 

The correlation functions of the atomic field, contrary 
to the electromagnetic waves, depend not only on the 
temperature, but also on the interactions and the dimen¬ 
sionality of the system. The density and phase of a three- 
dimensional condensate do not fluctuate. The situation 
is different in lower dimensions. In a two-dimensional 
Bose gas the Berezinskii-Kosterlitz-Thouless transition 
is manifested by a characteristic behavior of correlation 
functions m- In this paper we concentrate on a quasi- 
one-dimensional system, i.e. the situation when both: 
the chemical potential fi as well as the thermal energy 
ksT are smaller than a transverse confinement energy 
fi^ksT < hujr^ where ujr is the transverse trap frequency. 
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One-dimensional systems are different than their 2 D 
or 3D analogues. First of all there is no phase transition 
from the thermal to the condensate phase both in a trap 
and a uniform system. There exists however a degener¬ 
acy temperature j 2m (n is a ID density) below 

which the quantum effects become important in the uni¬ 
form gas at the thermodynamic limit. Similarly, for the 
ideal ID gas confined in the harmonic trap of frequency 
cj, the characteristic temperature of quantum degener¬ 
acy is equal to log(2A^). If, in addition 

to a harmonic confinement, the gas is weakly interact¬ 
ing, then there exist a second characteristic temperature 
= lb{hWz)‘^N/‘^2ii < Tc. As shown in [ 11 ] for temper¬ 
atures range in between the two <T <Tc^ the system 
is in a quasicondensate regime where density fluctuations 
are suppressed but phase fluctuations are large. Only 
below T(j) the system enters a fully coherent region with 
suppressed density and phase fluctuations. The phase 
fluctuations have been observed in expansion when they 
are transformed into density fluctuation m- The state¬ 
ment of suppressed density fluctuations remains a com¬ 
mon believe as there is no direct experimental evidence 
supporting or disclosing the above mentioned theoretical 
prediction based on the Bogoliubov approach. 

For a very elongated quasi-one-dimensional system the 
local density approximation is justified. Then both: the 
trapped and the uniform systems should exhibit simi¬ 
lar properties. The uniform system interacting via the 
zero-range pseudo-potential of the strength g is formally 
exactly soluble via the Bethe ansatz. In the limit of 
strong interaction, i.e. if 7 = mg/li?n ^ 1 the sys¬ 
tem reaches the Tonks-Girardeau limit of impenetrable 
bosons i.e. the regime of fermionization. The system we 
investigate here is in the opposite regime of parameters, 
namely the regime where interactions are small. 

One-dimensional Bose gases are special systems in a 
sense that, as discovered many years ago by Lieb m, 
they exhibit two families of elementary excitations. As 
already identified by Lieb na, the main branch is related 
to phonons. It turns out that dark solitons can be associ¬ 
ated to the type II branch of the elementary excitations. 
That has been demonstrated by analyzing the disper¬ 
sion relation for solitary waves m and also by studies 
of quantum nature of dark solitons beyond the mean- 
field [15]. Recently, we have shown that dark solitons are 
spontaneously generated in quasi-one-dimensional Bose 
gases at equilibrium m- By analyzing the statistical 
distributions of excitations within both the Lieb-Liniger 
model m and the classical fields approximation [ia[l9] 
we proved that type II excitations are indeed quantum 
solitons m- 

The spectrum of elementary excitations plays a cru¬ 
cial role for the correlation functions and coherence of 
the system. There is a reach literature devoted to both 
the first order ED HD and the second order correlations 
[23U26] of the ID or quasi-one-dimensional systems. A 
number of different approximate formulas for the second 
order correlation function in various regimes of the sys¬ 


tem parameters and temperatures are available. In par¬ 
ticular in the region of fermionization the antibunching 
of atoms is predicted. On contrary, at temperatures close 
to the degeneracy temperature and in the limit of weak 
interactions, the bunching of atoms is expected. The 
piece-wise valid formulas, although quite interesting, do 
not give a simple and clear picture of the density-density 
correlations in the entire range of temperatures across 
the ‘BEC-like transition’ for a weakly interacting system 
where the Bogloliubov theory adopted to elongated sys¬ 
tem can be compared to the predictions based on the 
Lieb-Liniger model. 

Observation of a spatial dependence of the second or¬ 
der correlation functions is difficult for several reasons. 
First of all observations in situ are limited by a finite res¬ 
olution of detectors and, in addition, the standard detec¬ 
tion techniques involve a partial spatial averaging - there¬ 
fore a column density is monitored only. The most com¬ 
mon method of detection is destructive and the atomic 
cloud is observed after releasing atoms from the trap fol¬ 
lowed by the expansion. All these effects smear out most 
of density fluctuations. Earlier studies of the second or¬ 
der correlation function in expansion considered ballistic 
expansion of noninteracting atomic cloud m showing 
that the correlation function observed is related to the 
correlation in situ by some scaling only. In [28] the au¬ 
thors showed that density-density correlations of the in¬ 
teracting uniform gas manifest themselves in the expan¬ 
sion in a form of density ripples. The spectrum of the 
ripples is related to density-density correlations. The re¬ 
cent experimental results of [29] for the two-point density 
correlation measured across the Bose-Einstein condensa¬ 
tion temperature in elongated expanding system are still 
awaiting full understanding and theoretical description. 
The authors support the experimental results with the 
ideal gas theory and evidently one should account for in¬ 
teractions to get a deeper insight into measured density 
correlations. 

In this paper we are going to demonstrate that the 
properties of the density correlations of a Bose gas at 
temperatures below the critical one, after the gas is re¬ 
leased from the trap, indicate large density fluctuations 
in the quasicondensate regime and might be a signa¬ 
ture of the presence of thermal solitons in the system. 
As our calculations show, for elongated Bose gases the 
density correlation function takes large values below the 
threshold temperature, i.e. for T < Tc. These values 
are larger than those calculated for a weakly-interacting 
Bose gas within the Bogoliubov approximation and dif¬ 
fer even more significantly from the results obtained for 
the ideal Bose gas theory. This indicates that not only 
phonons but also the second branch of elementary exci¬ 
tations influences the properties of density correlations. 
Large values of the density correlation function are there¬ 
fore a signature of the presence of dark solitons. This 
should not be very surprising as an every dark soliton 
is associated with a density dip. The recent experiment 
showing large values of the density correlation function 
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for temperatures below the critical one [53] (see Fig. 3c) 
support our finding. 

Below we are going to describe our numerical approach 
to determine the density-density correlation function of 
the elongated interacting Bose gas after expansion from 
the harmonic trap, taking into account a realistic de¬ 
tection procedure, in particular a column averaging and 
finite resolution of a CCD camera. 

In second quantization formalism one introduces the 
field operator ^(r) (^’^'(r)) which annihilates (creates) an 
atom at position r. The second order correlation func- 
tion, G(2)(r,r'), is then defined as a statistical average 


atomic clouds at the thermal equilibrium. While evolv¬ 
ing with the help of Eq. (§ the classical field reveals 
the properties of an atomic cloud within a single-shot ex¬ 
periment. We therefore pick up typically thousands of 
members of the canonical ensemble and propagate them 
without an external trapping potential, i.e. Vtrap = 0, 
for a given time. We then compute the autocorrelation 
function of the density integrated along one of the radial 
directions n 2 D{x^z) = / |T(x, z)p(i^, where 2 ; is the 

axial direction. However, because of the finite resolution 
of the optical imaging system we smooth the radially in¬ 
tegrated density via the following convolution: 


G(2)(r,r') = (4't(r)4r(r)4'1'(r')^(r')) • (1) 


It is convenient to use a normalized second order corre¬ 
lation function 


5(2)(r,r') = 


G(2)(r,r') 


G'(i)(r,r)G(i)(r',r') ’ 


( 2 ) 


n 2 _D(r) = ^ J ^ 2 d(u) exp [-(r - u)^/2cr^] dPu , (4) 

with the Gaussian function of the half width at 
half maximum (HWHM) of about a few micrometers 
and A = The density correlation func¬ 

tion itself is obtained by averaging the autocorrelation 
f n 2 D(u)n 2 n(u + r) d^u over all performed realizations: 


where G^^^(r, r') = (T’*'(r)T(r')) is the first order corre¬ 
lation function. 

To calculate the second order correlation function we 
turn to the approximate treatment of an interacting Bose 
gas, which is called the classical fields method (CFA) 
namj. The idea of the approach can be considered as 
an extension of the original Bogoliubov idea m, i.e. the 
bosonic field operator T(r) is replaced by the complex 
wave function T (r). This wave function corresponds to a 
state having large energy related to the temperature and 
is composed of macroscopically occupied single particle 
modes only. To obtain the thermal equilibrium state of 
an elongated weakly interacting Bose gas we generate 
an ensemble of classical fields T(r) corresponding to a 
given temperature. An effective way of getting members 
of such a canonical ensemble is to use the Monte Carlo 
algorithm |31| . 

Each classical field belonging to the canonical ensemble 
obeys the following equation of motion QSI: 










(3) 


where the coupling constant g describes the contact in¬ 
teraction between atoms. It looks like the usual Gross- 
Pitaevskii equation for the Bose-Einstein condensate at 
zero temperature. However, here the complex function 
T(r) carries information on both the condensed and non 
condensed atoms. The condensate and the thermal cloud 
can be split via the coarse-graining procedure mm- 
In short, the coarse-graining allows to separate the mode 
characterized by the largest first order coherence length 
(or a coherence time) from the remaining part of the clas¬ 
sical field. 

Different realizations of classical fields, being the mem¬ 
bers of the canonical ensemble, represent the set of 


G^'^\0,r) = {y n2D(u)n2D(u + r)d^u). (5) 

To get the normalized density correlation function we 
divide ^ b y the autocorrelation function of the mean 
density [33] . integrated along the radial direction: 

g^^)(0,r) = G^^)(0,r) / J{n 2 Diu)) {n 2 D{u + r)) (fu. 

(6) 

Autocorrelation functions appearing in both the numer¬ 
ator and denominator of the above definition are cal¬ 
culated using the fast Fourier transform. 

Because we use the classical fields approximation, i.e. 
we replace the bosonic fields operators by the complex 
functions, the shot-noise term does not appear in our ap¬ 
proach. The origin of the shot-noise term is in the com¬ 
mutation relations between the creation and annihilation 
operators for the bosonic field. Obviously in the classical 
fields approximation the bosonic commutation relations 
are not accounted for. It might be questionable to use 
the CFA for obtaining the density correlation function 
because the shot-noise term appears explicitly in the ex¬ 
pression for the second order correlation function while 
expressed with the help of the first order correlation func¬ 
tions. The shot-noise term is of the order of the number 
of particles in the system, N, whereas the other terms are 
proportional to N‘^. Hence, it can be neglected for large 
systems. Although, it is detectable in the experiment of 
Ref. [29|, a special care is taken to exclude the atomic 
shot-noise peak from the experimental data. 

In Fig. we show the density correlation function, 
z), for the temperature T = 3.15 nK after free ex¬ 
pansion over 46ms (upper and lower-left frames). The 
Bose gas consisting of 10^ rubidium atoms is initially 
confined in a quasi-one-dimensional harmonic trap with 
the radial frequency = 27r x 113 Hz and the axial one 
= 27r X 1 Hz. The aspect ratio is just as for the most 
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FIG. 1: (color online). Axial (solid lines) and radial (dashed 
lines) cuts of the density correlation function g^‘^\x,z) of 
the Bose gas after 46ms of expansion (upper frame). The 
Bose gas consisting of 10^ rubidium atoms is initially con¬ 
fined in a harmonic trap with radial and axial frequencies 
cjr = 27r X 113 Hz and = 27r x 1 Hz, respectively. The sys¬ 
tem is prepared in the thermal equilibrium at the temperature 
T = 3.15 nK. Three sets of curves are obtained by taking dif¬ 
ferent values of the smoothing parameter: cr = 4 /im, 6 gm, 
and 8 fim (from top to bottom). Lower panel: The density 
correlation function g^^\x,z) for a = 8 gm (left frame) and 
the two-dimensional density (i.e., density averaged along the 
direction of imaging) after expansion, first smoothed with the 
parameter cr = 8 gm and next averaged over all realizations 
(right frame). 


elongated trap used in the experiment of Ref. [29] . Sim¬ 
ilarly to the experiment, the theoretical system is also 
in the quasi-one-dimensional limit, i.e. its chemical po¬ 
tential is smaller than the radial excitation energy, hujr- 
Some other parameters are different, in particular the 
number of atoms is 10 times smaller. We use the system 
with smaller number of atoms confined in a trap with 
lower trap frequencies to decrease its mean field energy. 
In such a case the expansion of the gas after the trap 
is removed is much slower and becomes feasible numeri¬ 
cally. We believe, however, that the scaling does not alter 
significantly the main physical processes responsible for 
behavior of the density correlation function. 

The gas is then released from the trap and expanded 
for 46 ms. For that we solve the Eq. § with no trapping 
potential, Vtrap = 0- The expansion is not rapid and after 
46 ms the radial size of the cloud remains at the order of 
tens micrometers (see lower-right frame in Fig. where 
we show the density, first smoothed with the parameter 
a = 8 gm and then averaged over all realizations). As op¬ 


posed to the experiment of Ref. [HI in our calculations 
we image the whole radial extension of the cloud. The 
upper frame in Fig. presents the density correlation 
function cuts along the axial (solid curves) and radial 
(dashed curves) directions. The bunching at short inter¬ 
particle distance (g^^^ > 1) as well as the oscillations of 
correlations along axial direction with the values below 
one are clearly visible. Our results correspond to the ob¬ 
servations in the experiment of Ref. [Hj (see also earlier 
experimental data of Ref. [33]). At large distances which 
coincide with the edge of the atomic cloud, the density 
correlation function reaches its limit equal to one. 

Speaking more precisely, during the expansion, in the 
limit of large expansion time, we monitor the momentum 
distribution, i.e., the Fourier transform of the classical 
field. We find that it does not change in time almost at 
all what means that the interaction energy does not play 
any significant role during the expansion in our case, i.e., 
for the system with N = 1000 atoms. In practice the 
classical field evolves freely and at any time can be found 
with the help of the propagator of the free Schrodinger 
equation. This technique is much faster than direct solv¬ 
ing of the Eq. ^ on the grid and, as we checked, gives 
the same results. In fact, this technique was used by us 
since while calculating the density correlation function, 
the Eq. ^ has to be averaged over a few thousand of 
realizations. 

Next, we analyze the peak value of the second order 
correlation function, ^*^^^(0, 0), of expanding Bose gas for 
various temperatures. There is a huge discrepancy be¬ 
tween the experimentally measured values and the re¬ 
sults predicted by the ideal Bose gas theory for elongated 
systems, see Ref. [Hj- For temperatures below the criti¬ 
cal one the peak height of g^‘^\x^ z) rapidly decreases to 
unity for the ideal gas case whereas experiment proves 
that the peak height does not change significantly within 
a large range of temperatures. In Eig. we mark the 
results obtained within the classical field approximation 
by red dots. The temperature is given in units of the 
critical temperature for the ideal Bose gas corrected due 
to finite number of atoms and the quasi-one-dimensional 
character of the trapping potential. Eor the ideal Bose 
gas in the thermodynamic limit the critical temperature 
is Tc = hjujN^^^/ks [C(3)]^^^, where Cj = and 

CO is the Riemann zeta function [34] . Eor the parameters 
we use in our calculations, Tc = 10.5 nK. After correc¬ 
tions mentioned above the critical temperature is, how¬ 
ever, shifted down and equals Tc = 6.3 nK [34] . Since in 
the system we consider the number of atoms is very small, 
the interactions do not influence the critical temperature 
much. 

Now we compare our results with the ones obtained 
within the theoretical model described in Ref. [28] • This 
model allows to calculate the two-point density corre¬ 
lation function of very elongated ultracold Bose gas af¬ 
ter its release from the trap. Eor a tight radial confine¬ 
ment the transverse atomic motion is essentially frozen 
and thus decouples from the motion along the long axis. 
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FIG. 2: (color online). Upper frame: Peak height of the 
density correlation function as calculated within the classical 
field approximation (red dots), the ideal Bose gas model as in 
Ref. [27] (green solid line), and the Bogoliubov approximation 
extended to low-dimensional quasicondensates as in [28] (red 
dashed line). The critical temperature equals Tc = 6.3 nK and 
the smoothing parameter a = 8 /im. The shaded area shows 
the range of temperatures when deep spontaneous dark soli- 
tons appear in the system. Note large values of the density 
correlation function below the threshold temperature calcu¬ 
lated within the classical field approximation. Lower frame: 
ksTIji (blue dots) and K/tt (red squares) as functions of tem¬ 
perature to confront the validity of the model of Ref. [28] (red 
dashed line in the upper frame). 


clear that the interactions strongly influence the density 
correlation function. 

Fig. [2] clearly demonstrates the differences between 
theoretical predictions based on the ideal Bose gas model 
as well as on the model of an interacting gas within the 
Bogoliubov approximation and the ones obtained in the 
framework of the classical fields approximation. Only 
CFA results qualitatively recover the experimental data 
of Ref. [29] showing the significant density correlations 
below the critical temperature for the elongated Bose 
gases. The reason is that CFA reproduces the correct 
spectrum of elementary excitations including not only the 
Bogoliubov phonons but also the II type excitations - the 
dark solitons. They do occur spontaneously in a quasi- 
one-dimensional interacting Bose gas at equilibrium as 
it was already reported in m- The dark solitons do 
not affect the first order correlation function, i.e. the 
coherence length as shown in m, however they affect 
the density fluctuations. Large peak values of density 
correlations can be associated with the existence of soli¬ 
tons, see Fig. As it was already shown many years 
ago by Lieb muT] and later on generalized to the fi¬ 
nite temperature case by Yang and Yang [35], a one¬ 
dimensional system exhibits two families of excitations. 
The one branch belongs to phonons and is well described 
by the Bogoliubov approximation. This kind of excita¬ 
tions are included in the model of Ref. m- However, 
there is the second branch of excitations which are dark 
solitons m which is neglected in the considerations in 
Ref. [28]. Hence, the origin of the discrepancy between 
the two approaches is, in our opinion, related to the pres¬ 
ence of the second type of excitations in one-dimensional 
systems, i.e. to the dark solitons. They are especially 
important in the regime of a quasicondensate where very 
deep solitons are possible due to large thermal energy. 


Therefore, the problem of calculating the density correla¬ 
tions is effectively reduced to ID. The density correlation 
function turns out to be at any time analytically related 
to the spectrum of “density ripples” which, on the other 
hand, can be obtained within the Bogoliubov approxima¬ 
tion extended to low-dimensional quasicondensates. The 
shape of the density correlation function depends on the 
expansion time, temperature, and the ID density. The 
ID density is calculated as the average of the radially in¬ 
tegrated central densities of classical fields belonging to 
the canonical ensemble corresponding to the given tem¬ 
perature. The red dashed line in Fig. upper frame, 
is the output of such calculations. One has to remem¬ 
ber, however, that the validity of Ref. [28] ’s model is re¬ 
stricted to such temperatures that ksT/(i Kj it^ where 
K = irhniD / {me) is the Luttinger liquid parameter with 
c being the speed of sound. In Fig. lower frame, we 
show both quantities, ksT/jj. and K = tt, to visualize the 
above mentioned condition. Certainly, the results of this 
model are questionable for temperatures about the criti¬ 
cal one. On the other hand, for lower temperatures it is 
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FIG. 3: Solitons for trapped atomic clouds at thermal equi¬ 
librium for the temperatures T = 4.94 uK and T = 6.09 uK. 
The panels show the time evolution of the density along the 
axial direction. Dark solitons are clearly visible as dark lines 
winding from the one edge of atomic cloud to the other. 


Above the critical temperature the peak density corre¬ 
lation should be equal to two, indicating atom bunch¬ 
ing. The CFA calculations, however, give a value of 
the second order correlation function close to one rather 
then two. This is related to the fact that the character¬ 
istic length over which the particles bunching vanishes 
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(i.e. = 1) decreases with increasing temperature. 

Eventually, this length becomes smaller than the spatial 
resolution of the imaging system and the bunching ef¬ 
fect cannot be visible any more. Our calculations take 
into account the realistic resolution of detectors, there¬ 
fore they give a decrease of the correlations above the 
critical temperature as it is observed in the experiment 
[29]. In our case, the decrease of is slower since 

we have fewer number of atoms in the system. 

In summary, we have studied the density correlation 
function of expanding Bose gas. The weakly interact¬ 
ing Bose gas at thermal equilibrium is confined initially 
in a very elongated trap. The whole range of temper¬ 
atures is considered, from temperatures as low as those 
of a pure condensate through the ones typical for quasi¬ 
condensates up to temperatures much above the critical 
one. Below the critical temperature the normalized den¬ 
sity correlation function does not fall rapidly to the value 
one as expected for the coherent system and predicted 
for the ideal Bose gas m Our finding sheds a new 
light onto the understanding of coherence of elongated 
or quasi-one-dimensional system. Contrary to what was 
expected, relatively large phase and density fluctuations 
signify a lack of a full first and second order coherence 


and atom bunching in a large range of temperatures be¬ 
low the degeneracy temperature. Indeed, the interactions 
are responsible for these features. However, the theory 
based on the Bogoliubov approximation alone seems not 
to be able to fix the problem. This is because in elon¬ 
gated systems besides the phonon excitations described 
by the Bogoliubov approximation there exist another ex¬ 
citations as predicted by the Lieb-Liniger model. These 
are the dark solitons. Our approach captures both kinds 
of excitations and therefore qualitatively recovers the ex¬ 
perimental data of Ref. [29] . 
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